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Levy flights are random walks where the dynamics is dominated by rare events. Even though 
they have been studied in vastly diflFerent physical systems, their observation in a single quantum 
system has remained elusive. Here we analyze a periodically driven open central spin system and 
demonstrate theoretically that the dynamics of the spin environment exhibits Levy flights. For the 
particular realization in a single-electron charged quantum dot driven by periodic resonant laser 
pulses, we use Monte Carlo simulations to confirm that the long waiting times between successive 
nuclear spin flip events are governed by a power-law distribution; the corresponding exponent rj = 
— 3/2 can be directly measured in real-time by observing the waiting time distribution of successive 
photon emission events. Remarkably, the dominant intrinsic limitation of the scheme arising from 
nuclear quadrupole coupling can be minimized by adjusting the magnetic field or by implementing 
spin echo. 


Levy flights are continuous-time random walks that are 
characterized by diverging average waiting times or aver¬ 
age step lengths. Levy flights occur if either the waiting 
time At or the step length Ax distribution show a power- 
law dependence, P{At) ^ Af^ or P{Ax) ^ Ax'^, with an 
exponent —2 < 77 < — 1. The fat tail of the distribution 
leads to extreme events that dominate the dynamics [l|. 
Levy flights have been observed experimentally in vari¬ 
ous physical systems, such as light scattering in specially 
engineered materials [1] or in subrecoil laser cooling [3. 
In all experiments to date. Levy dynamics had to be in¬ 
ferred from ensemble properties; tracking Levy flights of 
a single quantum system in real time has so far remained 
elusive. 

In this Letter, we theoretically investigate an open cen¬ 
tral spin system driven by periodic laser pulses and show 
that the random walk of the polarization of the spin en¬ 
vironment can be characterized by Levy flights. To be 
specific, we derive the dynamics for a particular realiza¬ 
tion of the central spin model: an electron confined in 
a single self-assembled quantum dot (QD) that interacts 
with the nuclear spin environment. This choice is mo¬ 
tivated by striking experimental results on an ensemble 
of self-assembled QDs which are compatible with Levy 
flights of the Overhauser field [3^1| • The random walker 
in this central spin system corresponds to the nuclear 
Overhauser held Bqi,, the step length to the shift in the 
Overhauser held ABoh induced by a nuclear spin flip 
and the waiting time is given by the time interval Atnuc 
between two successive nuclear spin flips. Nuclear spin 
flips in this system are enabled by spontaneous emission 
from the laser excited central (electron) spin. Conversely, 
long waiting times are characterized by a corresponding 
reduction in the photon emission rate. We confirm that 
the waiting time has a power-law distribution with an 
exponent 77 = —3/2 leading to extremely long waiting 
times. These Levy flights can be observed by monitoring 
the waiting time distribution of successive photon detec¬ 
tion events. Owing to the favorable ratio of the rates of 
emitted photons to nuclear spin flips it is experimentally 


feasible to probe the nuclear spin dynamics in real-time 
even with a modest photon detection efficiency. 

We consider an electron in a single self-assembled QD 
coupled to the radiation field reservoir and illuminated by 
a (cr+) circularly polarized laser pulse train propagating 
along the QD growth direction x. The external magnetic 
field B = Bz is applied transverse to the light propaga¬ 
tion direction (Voigt geometry) as shown in FIG. [T] (a). 
The electronic ground states (|tx), |4.x)) and the excited, 
negatively-charged trion states (IRlIx), It-l'-D-x)) realize a 
four-level system with ground (e) and excited (t) states 
processing in the external magnetic field B with frequen¬ 
cies We = gehi-BB and Wt = g^g-BB, respectively. The ef¬ 
fective g-factors of the electron (hole) are given by ge < 0 
{gh > 0), and 713 is the Bohr magneton. We describe the 
system in the basis of the light propagation direction and 
hereafter suppress the subscript x on the basis states, i.e. 
Itx) = It)- The Hamiltonian is then given by {h = 1 ) 

i?o=+y(|t) (il + lt) (tl) 

+y(im) (tui + itu) (mi) 

+wte(im) (mi + im) (mi) 

where the last term corresponds to the optical transition 
energy wte. 

A cr+-polarized laser pulse couples the states |t) and 
Ittfl) 0 while leaving |t) unaffected. Hence, |t) is called 
bright and |t) dark state. The duration of the laser pulse 
tp is taken to be short enough so that the coherent evo¬ 
lution due to the four-level Hamiltonian Hq and sponta¬ 
neous emission are negligible during the pulse. We jus¬ 
tify later that all elements of the electron density matrix 
p related to the excited states are zero at the arrival of 
a pulse. Then the populations before (i) and after (f) a 
single pulse are given by [13 
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FIG. 1. (a) The magnetic field B is applied in z-direction and 
the laser pulse train propagates along the x-direction, which 
is the QD growth axis, (b) The electron and the trion precess 
in the external magnetic field (green) and the (T+-polarized 
laser only couples to I't'a:) (blue). Spontaneous emission is 
measured in the circularly polarized basis (gray). On one 
hand, the coupling to the nuclear spin environment modifies 
the electron ground state precession frequency, while on the 
other hand the slow nuclear spin dynamics is dominated by 
optically-assisted nuclear spin flips (red). 


where the pulse area 9 = 2 f fi(t) dt determines the ex¬ 
cited electron population and fi is the laser Rabi fre¬ 
quency. 

The combined coupling of the four-level system to a 
periodic laser pulse train and the radiation field reservoir 
asymptotically forces the electron into a periodic quasi¬ 
steady state. This occurs by sequential excitation of the 
bright state by a laser pulse, followed by spontaneous 
emission of a photon at a rate T and the accompanying 
decay into one of the two ground states. If the pulse sep¬ 
aration T satisfies r ^ the trion population decays 
completely within the pulse separation t, which justifies 
(HD. The asymptotic bright state population before 
the pulse is then given by M 



sin^(a;er/ 2 ) 
i -b sin^(weT/2) 


provided 9 = tt and We S> T. Similarly, the coarse-grained 


spontaneous emission rate is given by 
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If the bright state population at the arrival of the pulse 
vanishes, so does the optical absorption and hence the 
coarse-grained spontaneous emission rate Tspon- This oc¬ 
curs exactly if the electron precession frequency satisfies 
the synchronization condition 

27r^ 

We G —Z. ( 2 ) 

T 

The expansion of Tspon close to a synchronization condi¬ 
tion m with m G Z gives 



The fact that Tspon depends quadratically on We is essen¬ 
tial for the emergence of Levy flights. 

The coupling of the central spin to the environment 
is given by the Fermi-contact hyperfine interaction be¬ 
tween the QD electron and the nuclear spin environment 

i[ni,[il hm = ^{sJi+s+p_). 

The number of nuclear spins is iV r; 10^ — 10® and 
is the spin projection operator to the 2 -axis of the 
nucleus for 1 < j < N. The hyperfine coupling con¬ 
stant Aj depends on the overlap of the electron wave 
function with the j*® nucleus. Equal coupling of all 
nuclei to the electron leads to the homogeneous hyper¬ 
fine coupling constant A = ^ = where 

An is the hyperfine coupling constant of the material, 
assuming a single nuclear species. The first term of 
iZhf, which we denote by IToh, describes the Overhauser 
shift 0 ; the electron processes around the vector sum of 
the external field and the effective magnetic Overhauser 
field Boh = Ajilp/igePB), created by the nuclear 

spins. The associated Overhauser precession frequency is 
<^Oh = The second term of describes 

nuclear spin flip-flop processes that are suppressed at 
high magnetic fields due to the energy mismatch between 
the electron Zeeman energy (we) and the average nuclear 
Zeeman energy (wn), i.e. We oJn 0. 

In self-assembled QDs, large anisotropic strain gives 
rise to strong nuclear quadrupolar interactions where the 
underlying electric field gradients point primarily along 
the x-axis. This quadrupolar interaction in conjunction 
with the Overhauser field term constitutes the dominant 
mechanism for optically assisted nuclear spin flips in self- 
assembled QDs 0 ; we therefore neglect the flip-flop pro¬ 
cess from here on. However, we note that the flip-flop 
processes in IThf, an effective non-collinear interaction of 
the form A^^Szli 0,[H,[I3| or a hole-induced 
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non-collinear process [lB | would yield the same dynam¬ 
ics. The dominant nuclear spin dynamics is described by 
the Hamiltonians 


N 






N 


Hoh = ^ AjSzii, 


^ 2 P{p+l) 




i=i 
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where is the Zeeman splitting of the j**' nucleus, and 
Wq is the quadrupolar energy for 1 < j < iV. 

The effect of the electron-nuclei interaction on the 
coupled system is twofold. On the one hand, the elec¬ 
tron dynamics is modified by the Overhauser field which 
changes cOe to the effective electron precession frequency 
= We -l- woh- The dark state is then only reached at 
the arrival of a pulse if synchronization (l2|) is fulfilled for 
the effective precession frequency 
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(3) 


The modified electron dynamics results in a coarse¬ 
grained spontaneous emission rate Tspon that depends on 
w|® and thus on the nuclear Overhauser field Boh, and 
vanishes if the synchronization condition is satisfied. 

On the other hand, the nuclear dynamics is modified 
by optical excitation of the electron, which results in nu¬ 
clear spin flip assisted spontaneous emission and leads 
to a slow time evolution of the nuclear Overhauser field. 
To describe such processes, we apply a Schrieffer-Wolff 
transformation (l6l Il7j| to eliminate the non-energy con¬ 
serving terms of Hq. We find second-order nuclear spin 
flip processes where the additional energy is supplied by 
a spontaneously emitted photon (FIG. [T] (b)) occurring 
with a rate 0 


^^ 2 l sin^(wfT/2) 

T \+S\\P{ujfT/2) 






where e = a/TV under the assumptions of homoge¬ 
neous hyperfine and quadrupolar coupling, Wn ^ ^ and 
= ±5 for simplicity [l^. The rate of optically-assisted 
nuclear spin flips i?nuc is drastically reduced in the vicin¬ 
ity of a synchronization condition (j3|) . The signatures of 
this reduction can be observed in a single realization of a 
random walk of w|® {Iz)'- the random walk freezes near 
a synchronization condition which therefore constitutes 
a trap (red lines in FIG. [5] (a)). An unsynchronized elec¬ 
tron is optically excited and induces an unbiased random 



FIG. 2. (a) The random walk of the effective electron pre¬ 

cession frequency (blue) spends most of the time close to 
traps (red), indicating a Levy distribution for long waiting 
times. The inset resolves the random walk between the traps, 
(b) Tspon gives the spontaneous emission rate of the corre¬ 
sponding random walk in (a), (c) Exposing an ensemble of 
500 QDs with normally distributed initial nuclear polariza¬ 
tion {{Iz) = 5, (T = pN) (red) to the pulsed laser train for 
10®r = 10®/r leads to a final (blue) distribution of the effec¬ 
tive electron precession frequency that is heavily peaked at 
the different traps (inset). For all simulations we used the pa¬ 
rameters r = lO/F, We = 2802^, Wn = O.ISF, N = 2000, wq = 
O.OlSr, Ah = lOOF, and we assumed inhomogeneous hyper- 
fine coupling constants. 


walk of the nuclear polarization 1^ = Y^^=i{Pz) until w|® 
is synchronized. In the absence of other interactions, the 
time spent near such a trap can be on the order of the 
measurement time. We also simulated the time evolution 
of the distribution of w|® starting from a Gaussian distri¬ 
bution of (red histogram in FIG.[2](c)). We find 

that the probability of finding w|® within a single nuclear 
spin flip from a trap ([3]) exceeds 99%. These simulations 
are in excellent qualitative agreement with experiments 
carried out on an ensemble of single negatively-charged 
self-assembled QDs 

To quantify the freezing of a random walk of in the 
vicinity of a trap, we introduce the probability distribu¬ 
tion of the waiting time between two successive nuclear 
spin flip events Pnuc(Atnuc) Hi- Glearly, long waiting 
times occur for close to a trap since Atpuc ^ 

(wl® — ^m)“^. This relationship of Atpuc to a;|® allows 
us to write P„uc(Atnuc)dAtnuc = nnuc(w|®)da;|®, where 
the probability to reach an effective precession frequency 
wl® in the trap nnuc(w|®) is assumed to be uniform. Es¬ 
timating the change in Atnuc induced by an infinitesimal 
change in w|® via oc (w|® - ^m)-^ oc AtlLl, 

results in the power-law distribution 

f’nuc(At„uc) oc At~^P. (4) 

This expression is confirmed by a more rigorous 
derivation [I^ or 0 and is valid for At„uc ^ 

[13, [IS- The expectation value 

of this distribution does not exist. The Levy behavior of 
the nuclear spin flip waiting time is thus based on the 
quadratic dependence of Tspon and i?nuc on a;|®. 
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FIG. 3. The optical waiting time distribution for Ro — OF 
(red) is extracted from 100 time traces. The corresponding 
power-law exponent in this simulation is —1.49, consistent 
with the theoretical prediction of —3/2. For — 10®/F 
(blue) we find a regime where a power-law with exponent 
— 1.53 is observed, followed by an exponential decay. All other 
parameters are the same as in FIG. [2] 

The physical observable in an experiment is the wait¬ 
ing time Atgpon between two successive photon emission 
events, irrespective of whether the photon emission is ac¬ 
companied by a nuclear spin flip or not. The rate of 
photon emission events is given by Tspon and an anal¬ 
ogous derivation yields the waitin g ti me distribution of 
successive photon emission events [l0| 

^spon(Atspon) OC Algp^^j , (5) 

for Atspon » (^nucLeft^ 2 ^^±^) . As a conse¬ 

quence of the proportionality of the rate of nuclear spin 
flip events i?nuc to spontaneous emission events Fspon, 
the power-law exponents of the two waiting time distri¬ 
butions Pnuc and Pspon are identical. Since the rate of 
photon detection events is enhanced by a factor com¬ 
pared to the nuclear spin flip rate, it is possible to moni¬ 
tor the Levy flight behavior of nuclear trapping times in 
real-time even with modest photon detection efficiencies. 

To numerically verify the power-law behavior of the 
waiting time distribution, we averaged over 100 Monte- 
Carlo simulations of the type shown in FIG. [2] (a) and 
extracted a power-law exponent ry = —1.49 (red dotted 
line in FIG. [3|) which is in excellent agreement with the 
predicted distribution of optical waiting times ([5]). 

Experimentally observable waiting times are limited 
due to additional interactions that have been neglected 
so far. To investigate their impact on the probability 
of nuclear trapping times, we assume that they can be 
modeled as an additional constant nuclear spin flip rate 
Rq such that the overall nuclear spin flip rate becomes 

Pnuc — Pnuc + Ro- 

The constant rate Pq prevents diverging average trapping 
times since it effectively limits the trapping times to P^ . 
This scenario corresponds to a truncated Levy flight . 
We derived the modified distribution of nuclear trapping 
times [inl | 

Pnuc(At„uc) OC At-3/2(l -b 2i?oAt„uc)e“^‘>^‘““‘=. 


If Pq = 0 then expression o is recovered, but if 
Atnuc ^ Rq^ the exponential term dominates and leads 
to an exponential decay. However, in an intermediate 
regime, where the background rate Pq is much smaller 
than Pnuc at a single spin flip away from synchroniza¬ 
tion, i.e. Pq < PnucLetf=m 22 L+ 4 , a power-law regime 
is observed before the exponential cutoff. The optical 
waiting time distribution Pspon with constant offset Pq 
is again closely related to the nuclear trapping time dis¬ 
tribution Pnuc- A model simulation including a constant 
offset Pq ^ OF shows a power law distribution with coef¬ 
ficient — 1.53 (blue dotted line in FIG.[3l) for intermediate 
waiting times and an exponential decay for large waiting 
times. 

The physical origin of Pq can be classified as system- 
specific or intrinsic to the scheme. In the case of an 
electron localized in a self-assembled QD the coupling of 
the electron spin to phonons or to a Fermionic reser¬ 
voir [ 2 ^ could limit the trapping times. In addition, fi¬ 
nite pulse duration and imperfect laser polarization lead 
to a finite decay rate out of the electronic dark state. 
These effects can be suppressed by proper engineering. 
In contrast, the scheme can be intrinsically limited by 
the interaction that causes the optically assisted random 
walk, which is the quadrupolar interaction Hq for self- 
assembled QDs. In the time windows that separate the 
successive laser pulses, Hq leads to a coherent precession 
between different R eigenstates; the precession frequency 
is dominated by the energy detuning Wn- In the ab¬ 
sence of a measurement interaction interrupting the co¬ 
herent evolution no actual nuclear spin flip can occur due 
to the energy mismatch of the coherently coupled nuclear 
states. However, the laser pulse terminating the period of 
coherent evolution under Hq effects such a measurement 
and could lead to a nuclear spin flip. Gonsequently, the 
rate at which these intrinsic nuclear spin flips occur de¬ 
pends on the bright population in a different nuclear spin 
manifold at the arrival of a pulse. Similar intrinsically 
limiting mechanisms appear for systems with dominant 
non-collinear or hyperfine flip-flop interactions, where the 
latter is discussed in Q. 

The intrinsic contribution to Rq can be minimized by 
tuning the magnetic field B and pulse repetition period t. 
The coherent evolution of the nuclear spins is dominated 
by the nuclear Zeeman frequency cUn and can be reduced 
if the nuclei are also synchronized [l^ 

27r 

Wn e —z. (6) 

r 

In a system with several nuclear species it could be chal¬ 
lenging to synchronize all nuclei for an experimentally 
achievable B and r. Adding an echo pulse on the elec¬ 
tron would add flexibility to synchronize the nuclei and 
the electron separately. In this case, it is important to 
displace the echo pulse by a finite Afgcho from r/2 to 
ensure that the electron dynamics responsible for the op¬ 
tically assisted nuclear spin flips at a rate i?nuc is not 
completely reversed. The displaced echo pulses define a 
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new synchronization condition for the electron which is 
given by 


^techo 

This procedure adds an additional parameter Atgcho that 
controls the electron dynamics, while r can be extended 
to accommodate better synchronization for multiple nu¬ 
clear species without affecting the electron dynamics. 

In summary, we have demonstrated that the polariza¬ 
tion of the spin environment of a central spin system that 


is driven by a resonant periodic pulse train exhibits Levy 
flights. Thanks to the decoupling of the central spin and 
the environment and the proportionality of the environ¬ 
mental spin flip and spontaneous emission rates, our cal¬ 
culations indicate that the signatures of Levy flights can 
be observed by measuring the waiting time distribution 
of the emitted photons. While we illustrated the scheme 
with self-assembled quantum dots, we expect our findings 
to be relevant for a broad class of solid-state emitters in¬ 
cluding NV centers. More generally, our results show 
that the previously established connection between Levy 
flights and dark states S[i3 can be extended to systems 
with explicit time-dependence. 


[1] M. F. Slilesinger, G. M. Zaslavsky, and J. Klafter, 
Nature 363, 31 (1993) 

[2] P. Barthelemy, J. Bertolotti, and D. S. Wiersma, 
Nature 453, 495 (2008) 

[3] F. Bardou, J. Bouchaud, O. Emile, 

A. Aspect, and C. Cohen-Tannoudji, 

Physical Review Letters 72, 203 (1994) 

[4] A. Greilich, R. Oulton, E. Zhukov, I. Yugova, 

D. Yakovlev, M. Bayer, A. Shabaev, A. Efros, 

I. Merkulov, V. Stavarache, D. Reuter, and A. Wieck, 
Physical Review Letters 96, 227401 (2006) 

[5] a. Greilich, A. Shabaev, D. R. Yakovlev, A. L. Efros, 

I. a. Yugova, D. Reuter, a. D. Wieck, and M. Bayer, 

Science (New York, N.Y.) 317, 1896 (2007) 

[6] E. Barnes and S. E. Economou, 

Physical Review Letters 107, 047601 (2011) 

[7] S. E. Economou and E. Barnes, 
Physical Review B 89, 165301 (2014) 

[8] 1. A. Yugova, M. M. Glazov, E. L. Ivchenko, and A. L. 
Efros, Physical Review B 80, 104436 (2009), 0905.1254 

[9] B. Urbaszek, X. Marie, T. Amand, O. Krebs, 
P. Voisin, P. Maletinsky, A. Hogele, and A. Imamoglu, 
Reviews of Modern Physics 85, 79 (2013) 

[10] M. Issler, J. Holler, and A. Imamoglu, Supplementary 
material in preparation. 

[11] 1. T. Vink, K. C. Nowack, F. H. L. Koppens, 

J. Danon, Y. V. Nazarov, and L. M. K. Vandersypen, 
Nature Physics 5, 764 (2009) 

[12] C. Latta, A. Hogele, Y. Zhao, A. N. Vamivakas, 
P. Maletinsky, M. Kroner, J. Dreiser, 1. Garusotto, 


[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 


[ 21 ] 

[ 22 ] 

[23] 


A. Badolato, D. Schuh, W. Wegscheider, M. Atatiire, 
and A. Imamoglu, Nature Physics 5, 758 (2009) 

A. Hogele, M. Kroner, C. Latta, M. Glaassen, 

1. Garusotto, G. Bulutay, and A. Imamoglu, 
Physical Review Letters 108, 197403 (2012) 

G. Bulutay, Physical Review B 85, 115313 (2012) 

W. Yang and L. J. Sham, 

Physical Review B 85, 235319 (2012) 

1. Shavitt and L. T. Redmon, 
The Journal of Ghemical Physics 73, 5711 (1980) 

M. Issler, E. Kessler, G. Giedke, S. Yelin, 

1. Girac, M. Lukin, and A. Imamoglu, 

Physical Review Letters 82, 3008 (2010), 1008.3507 
J?nuc is derived for spin-1 particles by replacing —>■ 

this captures the relevant nuclear dynamics. 

F. Bardou, J. Bouchaud, A. Aspect, and G. Gohen- 
Tannoudji, Levy Statistics and Laser Cooling (Gam- 
bridge University Press, 2003). 

For nnuc(t^S^) to be uniform, the distance to the syn¬ 
chronization condition should not exceed a single nuclear 
spin flip A, where the waiting time is given by ri. There¬ 
fore, O only describes the tail of the distribution, i.e. 


At 

nuc S> n. 

R. Mantegna and H. 



T. Muller, A. Badolato, and A. 
Physical Review B 77, 075317 (2008) 


Stanley, 

1. Fal’ko, 

Galland, 

Imamoglu, 





